MULTI- VARIABLE POLYNOMIAL SOLUTIONS TO PELL'S 
EQUATION AND FUNDAMENTAL UNITS IN REAL 
§ ! QUADRATIC FIELDS 
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Abstract. Solving Pell's equation is of relevance in finding fundamen- 
tal units in real quadratic fields and for this reason polynomial solutions 
(f) , are of interest in that they can supply the fundamental units in infinite 

families of such fields. 

In this paper an algorithm is described which allows one to construct, 
for each positive integer n, a finite collection, {Fi}, of multi- variable 
■ polynomials (with integral coefficients), each satisfying a multi- variable 

polynomial Pell's equation 



Cf -FiHf = (-l) n -\ 



where d and Hi are multi-variable polynomials with integral coeffi- 
cients. Each positive integer whose square-root has a regular continued 
fraction expansion with period n + 1 lies in the range of one of these 
polynomials. Moreover, the continued fraction expansion of these poly- 
nomials is given explicitly as is the fundamental solution to the above 
■ multi-variable polynomial Pell's equation. 

Some implications for determining the fundamental unit in a wide 
class of real quadratic fields is considered. 



o 
o 
o 

^ ■ 1. Introduction 

-i— > 

Solving Pell's equation is of relevance in finding fundamental units in real 
quadratic fields and for this reason polynomial solutions are interesting in 
that they can supply the fundamental units in infinite families of such fields. 
There have been several papers written over the past thirty years which 
jjj ■ describe certain polynomials whose square roots have periodic continued 

fraction expansions which can be written down explicitly in terms of the 
coefficients and variables of the polynomials. See for example the papers of 
Bernstein Q, Levesque and Rhin |4|, Madden ||, Van der Poorten [1C] and 



Van der Poorten and Williams |11]. 



In this paper an algorithm is described which allows one to construct, for 
each positive integer n, a finite collection of multi-variable Fermat-Pell poly- 
nomials which have all positive integers whose square-roots have a continued 
fraction expansion of period n+1 in their range. If F{ := Fi(to, t±, ■ ■ ■ , t^n+i j ) 
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is any one of these polynomials, the fundamental polynomial solution to the 
equation 

(1.1) C f - F t Hf = i-ir- 1 

(where Cj and Hi are polynomials in the variables to, ti, ■ ■ ■ , t p+i j) can be 

found. Moreover, the continued fraction expansion of \fFl can be written 
down when t%, ■ ■ ■ , t, n +i i > and to > 9i{pli • ■ ■ , tin+ii), a certain rational 
function of these variables. Some implications for single-variable Fermat- 
Pell polynomials are discussed as are the implications for writing down the 
fundamental units in a wide class of real quadratic number fields. 

Definition: a multi-variable polynomial 

F:=F(t ,t lr -- ,t k ) eZ[t ,t u --- ,t k ], k> 1 
is called a multi-variable Fermat-Pell polynomial [] if there exists polynomials 

C := C(to, h, ■ ■ • , tfc) and H := H(tQ,tx, ■ ■ ■ ,t k ) G Z[tQ, ti, ■ ■ ■ , % ] 
such that either 

(1.2) C 2 — F H 2 = 1, for all ti, < i < k, or 

C 2 - FH 2 = -1, for all t h < i < k. 

Such a triple of polynomials {C, H, F} satisfying equation ( ^3 ) constitute 
a multi-variable polynomial solution to Pell's equation. 

Definition: The multi-variable Fermat-Pell polynomial F (as above) 
is said to have a multi-variable polynomial continued fraction expansion 
if there exists a positive integer n, a real constant T, a rational func- 
tion g(ti, ■ ■ ■ ,tk) G Q(ii, • • • , £fc) and polynomials ao := ao(to ; *i> • " > *fc) G 
Z[i Qj ti, • • • andaj := Oj(ti, • • • ,t k ) G Z[*i, ••• ,t k ], 1 < j < n, which 
take only positive integral values for integral ti > T,l < i < k and (possi- 
bly half-) integral to > 9(ti, ••• ,t k ) such that 

V~F = [oo; ai, • • • j a n, 2ao], for all ti 's in the ranges stated, < i < k. 
Remarks: 

(1) From the point of view of simplicity it would be desirable to replace 
the condition to > g(t\, ■ ■ ■ , tk) by to > T but it will be seen that for the 
polynomials examined here that the former condition is more natural and 
indeed cannot be replaced by the latter condition. 

(2) The restriction that the aj(ti,--- , t k ) > 0, 1 < i < n may also seem 
artificial to some since negative terms can easily be removed from a continued 
fraction expansion (see, for example [pX| D but this changes the period of the 
continued fraction so is avoided here. 

(3) It may also seem artificial to have ao depend on a variable to while the 

1 These polynomials are called "Fermat-Pell polynomials" here to avoid confusion with 
"Pell Polynomials" and also because Fermat investigated the "Pell" equation. 
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other Oj's do not but this will also be seen to occur naturally. 

(4) Finally, allowing to to take half-integral values in some circumstances 

may also seem strange but this also will be seen to be natural and indeed 

necessary. 

Definition: If, for all sets of integers {t' ,t'i,--- ,t' k } satisfying t' > 
g(t[,--- ,t' k ) and t\ > T, 1 < i < k, 

X = Ci(t' ,ti, ■ ■ ■ , t' k ), Y = Hi(t' ,t 1 , • • • , t' k ) 

constitutes the fundamental solution (in integers) to 

x 2 -^,*;,... ,t' k )Y 2 = (-i)^ 1 

then (Ci(to, t±, ■ ■ ■ , t^), Hi(to, t±, ■ ■ ■ , t^)) is termed the fundamental polyno- 
mial solution to equation Ql.l|). 



Standard notations are used: 

111 1 

o H ■ ■ ... — := a + 

ai+ a2+ 03+ on 



a\ + 



1 

a 2 H 



1 

a 3 _| 

a N 

To save space this continued fraction is usually written [oq; a\ , • • • , a n ] . The 
infinite periodic continued fraction with initial non-periodic part a$ and 

periodic part oi, a n , 2oo is denoted by [ao;ai, ,a n ,2ao]. The i-th 

approximant of the continued fraction [ao; a%, ■ ■ ■ , ] is denoted by Pi/Qi- 

Repeated use will be made of some basic facts about continued fractions, 
such as: 

(1-3) P n Q n -i — P n -iQ n = (— l) n 1 , 

Pn+l = O-n+lPn + Pn—li 
Qn+1 = a n+lQn + Qn-1, 

each of these relations being valid for n = 1, 2, 3 • • • . 

Before coming to the main problem, it is necessary to first solve a related 
problem on symmetric strings of positive integers. 

2. A PROBLEM CONCERNING SYMMETRIC SEQUENCES 

Question: For which symmetric sequences of positive integers ai,...,a n 
do there exist positive integers ao and D such that 

(2.1) v A D = [o ;oi, ,a n ,2a ]? 

Let Pi I Qi denote the ith approximant of the continued fraction 
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By the well known correspondence between convergents and matrices 

'0 1\ fax 1\ fan l\ (P n Pn-l 

1 O il VI 0) \Q n Q n -i, 







h- 1 





(; 










1\ fax 1\ a n 1\ /0 1\ /Pn-l Pn 

1 O il 1 1 \ 1 J 1 Q„_i Q„ 



Since the left side in the second equation is a symmetric sequence of sym- 
metric matrices it follows that 



(2.3) P n = Q n _ 



Suppose y/ D = [ao; a±, , a n , 2oq] = ao+/3, where (3 =[0; ax, , a n , 2ao] 

so that 

= [0;oi, ,an,2a + /3], 

(2a + /9)P„ + Pn-i /3P„ + (2a P„ + P„-i) 



/3 



(2a + /9)Q„ + Qn-l PQn + (2a Qn + Qn- 



P 2 Q n + (2a Qn + Qn-l " P„)/3 - (2a Pn + Pn-l) = 0, 



(3 2 Q n + (2a Qn)/3 - (2a Pn + P n -i) = 0, (by (O)) 



/r> i/3 / 2 , 2aoP n + Pi-1 

The problem now becomes one of determining for which symmetric se- 
quences of positive integers ai,...,o n does there exist positive integers ao 
such that (2aoP n + P n -x)/Qn is an integer. 

Theorem 1. There exists a positive integer ao such that (2aoP n +P n -i)/Qn 
is an integer if and only if P n -\Q n -\ is even. 



Proof. Suppose first of all that P n _xQn-x is even. By equation ( |1.3| ) 

PnQn-l + ( — 1)™ = Pn-xQn 

(i) Suppose n is even. Then P n Q n ^\P n ~\ + Pn-i = P%-xQn- Choose t 
to be any integer or half-integer such that tQ n is an integer and ao := 
Qn-lPn-l/2 + tQ n > 0. Then 

2«pPn + Pn-l _ Qn-lPn-xPn + 2tP n Q n + Pn-l _ r^p _j_ p2 
Qn Qn 

(ii) Similarly, in the case n is odd, -P n Q n -\P n -x + Pn-i = -P^_i<5n- 
Choose t to be any integer or half-integer such that tQ n is an integer and 
a o := — Qn-iPn-i/2 + tQ n > 0. In this case 

2q P n + P n _i _ n+r> 2 

^Sn 

=>- Suppose next that P n _i and Q n -\ are both odd and that there exists a 
positive integer ao such that (2a P n + Pn-i)/Qn is a positive integer, m, say. 
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n 



Using ( |1,3| ) and ( p.3[ ) it follows that Q n is even. Then 2aoP n + P n -\ = mQ 
implies P n -i is even - a contradiction. □ 



Remarks: 

(i)Note that this process gives all ao such that (2aoP n + P n __i)/Q ra is an 
integer. Indeed, 

(2a P n + P n -i)/Qn = k, an integer 

<^=^ 2aoP n Q n -i = —P n -\Q n -\ + kQ n Q n -i, 

^ 2a (-l) n - 1 = 2a (P„Q„_i - P n -iQ n ), ( by Q) 

= —Pn-lQn-l + Qn(kQ n -l ~ 2aoP n -l), 

, ,-, n -l { —Pn-lQn-l . ^ kQ n -l — 2aoP n -l 
^=>- «0 = { — tj ^ hVn 



Notice also that if there is one such ao that there are infinitely many of 
them. 

(ii)Notice that, with P n , P n _i,Q n and Q n -i as defined above, if there 
exists a positive integer D satisfying (|2.l|) then D = p(to), for some allowed 
to, where 

P(t) = ( 9 "- 1 2 P - 1 + tQ B ) 2 + 2tP n + P^, t > ~ Qn 2 Q n n - 1 , (neven), 

'Qn— lPfi— 1 i \ mn 2 j. ^ Qn— 1 



The above theorem suggests a simple algorithm for deciding if, for a given 
symmetric sequence of positive integers ai, ■ ■ ■ , a n , there exist positive in- 
tegers ao and D such that (2.1) holds. Notice that all that matters is the 
parity of the en so all calculations can be done in Z2. First of all define the 
following matrices: 

J= (l 0)' *=(l J) and/ =(j 1 

Convert the sequence a\,a,2,--- ,a n to a sequence of J- and ET-matrices, 
according to whether each a$ is odd (replace by a X) or even (replace by a J). 
Prefix a J-matrix (to account for the initial in the continued fraction fl2.2| ) ) . 
Multiply this sequence together (modulo 2) using the facts that J 2 = K A = 
I, and JiT = if 2 J. 

/* A 

The final matrix = | J mod 2 there do not exist positive inte- 



gers ao and D such that (2.1) holds. 

Example 1. Do there exist positive integers ao and D such that 
VB = [a ; 22, 34, 97, 32, 15, 17, 17, 15, 32, 97, 34, 22, 2a ]? 
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As described above convert the sequence 22, 34, 97, 32, 15, 17, 17, 15, 32, 97, 
34, 22 to a sequence of J- and if-matrices, prefix a J-matrix and multiply 
the sequence together: 

J J JKJKKKKJK J J = JK(JK)JK = J K(K 2 J) J K 
Therefore there do exist positive integers ao and D such that 



D = [ao; 22, 34, 97, 32, 15, 17, 17, 15, 32, 97, 34, 22, 2a ]. 

3. Multi-variable Fermat-Pell Polynomials 
Definition: If { a±, ■ ■ ■ , o n } is a symmetric zero-one sequence such that 

(!;)n(? ;K — 

then the sequence { 01, ■ • • , a n } is termed a permissible sequence. Let r(n) 
denote the number of permissible sequences of length n. 

Note: It is not difficult to show that r(2m) = ((-l) m + 2 m+1 )/3 and that 
r(2m + 1) = ((-l) m + 5 x 2 m )/3. 



If D is a positive integer such that yD = [ao;a\, ,a n ,2ao] then 

{ oi, • • • , a n } mod 2 must equal one of the above permissible sequences and 
D is said to be associated with this permissible sequence . The collection 
of all positive integers associated with a particular permissible sequence is 
termed the parity class of this permissible sequence. Sometimes, if there is 
no danger of ambiguity, these collections of positive integers will be referred 
to simply as parity classes. 

Theorem 2. (i)For each positive integer n there exists a finite collection of 
multi-variable Fermat-Pell polynomials Fj(to, t%, ■ ■ ■ , ti n+ii), 1 < j < r{n), 
such that each positive integer whose square root has a continued fraction 
expansion with period n + 1 lies in the range of exactly one of these polyno- 
mials. Moreover, these polynomials can be constructed; 

(ii) These polynomials have a polynomial continued fraction expansion 
which can be explicitly determined; 

(Hi) The fundamental polynomial solution 
C = Cj(to,ti, ■ ■ ■ ,^n±ij), H = Hj(to,t\, ■ ■ ■ , t^n+ij) to 

(3.1) CP-Fjfah,- - ,\^)H 2 = {-i^- 1 

exists and can be explicitly determined. 

Proof. (i)The proof will be by construction. 

Step 1: Find all permissible sequences. This will involve checking 2'- _ 2 _ J 
zero-one sequences in a way similar to the example (|l]) above. 
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Step 2: For each permissible sequence {a\,--- ,a n } create a new sym- 
metric polynomial sequence {a\(ti), 02(^2), •• • ,a n -i(t2),a n (ti)} by replac- 
ing each ai and its partner a n+ \-i in the symmetric sequence by aj(ij) = 
a n +i-i(U) = 2tj + 1 if ^ = 1 and by Ojfo) = a n+ i_j(ij) = 2tj + 2 if 
a, = 0. This new sequence will sometimes be referred to as the sequence 
{a\, ■ ■ ■ ,a n }, if there is no danger of ambiguity. Each of the integer vari- 
ables ti (in the polynomial being constructed) will be allowed to vary inde- 
pendently over the range < t{ < 00 and each of the new a^s will keep the 
same parity and stay positive. 

Step 3 As in ([T^), form the continued fraction 

11 11 

oi(ti)+ a 2 (i 2 )+' ' a n -i{t 2 )+ a n (ti) 
and calculate P n , Q n , P n -i and Q n -i for this polynomial continued fraction, 
where these expressions are now polynomials in the t^s. 

Step 4 Construct Fj := Fj(to,ti, ■ ■ ■ , iin+ii), the multi-variable Fermat- 
Pell polynomial corresponding to the particular parity sequence under con- 
sideration. This is simply done by defining 



(3.2) Fj 




""f"- 1 + t Q n ) 2 + 2t P n + Pl_ x , (n even) 



— +W» + 2t P n -P%_ 1 , (nodd) 



where (— l) n+1 Q n -iP n -i/(2Q n ) < to < 00 and to can take half-integral 
values if Q n is even and otherwise takes integral values. 

Every positive integer whose square root has a continued fraction expan- 
sion with period n + 1 lies in the range of exactly one of these polynomials. 
That these polynomials are multi-variable Fermat-Pell polynomials follows 
from equation ( |3.4| ) below. 

(ii) With to in the range given, then 



s/Fj — [ao(to,ti, ■ ■ ■ , ip+ij); ai{ti), ■ ■ ■ ,a n (ti), 2ao(to,ti, ■ ■ ■ , t^n+ij)], 
for all ti > 0. Here 



-iP„-i 



2 

-1 
2 

(iii) Notice (using (O) and (O)) that 



+ toQ n , (n even) 



(3.3) a = a (to,h,- ■ ■ ,W,) := < - Qn _ lPn _ 1 ' , . 

2 J 1 n 1 + t Q ra , (nodd). 



(3.4) (a Qn + Pn) 2 - (ag + (2a P n + P n ^)/Q n )Q 2 n = (-l)"" 1 . 

To see that (aoQn + -fn> Qn) is the fundamental solution to (|3.1| ), notice that 



= [<*o(*o> *li • • • i • • • , 2ao(to, *i, ■ ■ ■ , t^n+l j )]■ 

This has period n + 1 and the nth approximant is ao + P n /Qn = («oQn + 
Pn)/Qn and by the theory of the Pell equation (aoQn + Pm Qn) is the fun- 



damental solution to (3.1). □ 
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As regards fundamental units in quadratic fields there is the following 
theorem on page 119 of 0: 

Theorem 3. Let D be a square-free, positive rational integer and let K = 
Q(\/D). Denote by ep the fundamental unit of K which exceeds unity, by 
s the period of the continued fraction expansion for \[D , and by P/ Q the 
(s — 1 )-th approximant of it. 

If D ^ 1 mod 4 or D = 1 mod 8, then 

e = P + Q\Tb~. 

However, if D = 5 mod 8 ; then 

e = P + QVd. 

or 

e 3 = P + Q\[~D. 

Finally, the norm of eo is positive if the period s is even and negative oth- 
erwise. 

It is easy, working modulo 4, to determine simple conditions (on to ) which 
make Fj = 2 or 3 mod 4 and thus to say further, for a particular set of 
choices of t%, ■ ■ ■ , i^n+i j and for all odd or even to, that if Fj is square-free, 

then aoQ n + Pn + \/F~jQn is the fundamental unit in Q[ y/Fj ]. For example, 
suppose that n is even and that the original Q n -i determined from the 
permissible zero-one sequence is also even (so that P n -\ and Q n are both odd 
and P n = Qn-i is even). Then the multi- variable form of Q n ~i evaluated 
in Step 3 will also have all even coefficients. Suppose - n 2 ~ 1 = cq + ^ tj 
mod 2. (Here Co may be and the sum ]P t ■/ may contain some, all or none 
of the tj's ) It is easy to see that Fj(t , t\, ■ • ■ , i^n+i , ) = (co + Yl V +^o) 2 + 1 
mod 4. Even more simply, if the original original Q n -\ as in Step 1 is odd 
(here also the case n is even is considered) then P n -\ as evaluated in Step 
3 is even and it is not difficult to show that in fact P ra _i = 2 mod 4 (since 
for n even P n Q n _i — P n _\Q n = —1) and that Q n is odd, which leads to 
Fj(to, h, • • • , ti ra+ij) = ip + 1 mod 4. Similar relations hold in the case 
where n is odd. 

The polynomials constructed in theorem (^) take values in only one parity 
class, if all the variables are positive. However, given any two parity classes, 
there are multi-variable Fermat-Pell polynomials that take values in those 
two classes. 

Theorem 4. Let n be any fixed positive integer large enough so that the set 
of positive integers whose square roots have a continued fraction expansion 
of period n + 1 can be divided into more than one parity class. 

(i) Given any two parity classes of integers whose square roots have con- 
tinued fraction expansions of period n + 1, there are multi-variable Fermat - 
-Pell polynomials, which can be constructed, that take values in both parity 
classes; 
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(ii) These polynomials have a polynomial continued fraction expansion 
which can be explicitly determined; 

(Hi) If F = F(to,c,ti, ■ ■ ■ , t, n+i j) is any such polynomial then the funda- 
mental polynomial solution 

C = C(t , C, t\, ■ ■ ■ , t^n±l^),H = H(t , C, t\, ■ ■ ■ ,£|2±lj) 

to 

(3.5) C 2 - FH 2 = (-I)™" 1 
can be explicitly determined. 

Proof. As in Step 2 in theorem (§) a polynomial sequence {a\,--- ,a n } is 
created. Suppose L\ = {b\ ■ ■ ■ ,b n } and L2 = {ci,--- ,c n } are the per- 
missible sequences associated with the two parity classes. Let i%, ■ ■ ■ , i& be 
those positions < [^^-\ at which the sequences agree. For each of these 
v's set a ir (t ir ) = a n+ i-i r (t ir ) = 2t ir + 1, if c ir is odd and set a ir (t ir ) = 
o-n+i-iriUr) = 2tj r + 2, if Ci r is even. Subdivide the remaining positions 
(those positions < L^^J a ^ which L\ and L2 differ) into two subsets: those 
at which L\ has a and L2 has a 1 and those at which L\ has a 1 and L2 
has a 0. 

Suppose ij is a position of the first kind. Let ^.(c,^) = a n+ i_j (c, tj ) 
= c + 2 + 2ij.. Repeat this for all the positions ij in this first set. Likewise, 
Suppose ij is a position of the second kind. In this case let Gtj.(c, ij.) = 
q„ + i -ijfc, ii n+ij) = c + 1 + 2t^.This is also repeated for all the positions ij 
in this second set. Step 3 and Step 4 are then carried out as above. The rest 
of the proof is identical to theorem (Q) . Denote the polynomial produced by 

(3.6) F := F(to, c, ti, • • • , ^n±lj)- 

As in theorem (|2j) , if c and all the £j's are non-negative, 1 < i < L^^^J an d 
t > (-l)"+ 1 Q n „ 1 P n _ 1 /(2Q n ) then 

v A F= [ao;ai, ,a n ,2a ], 

where the a^'s, 1 < i < n are as defined just above and oq is as defined in 
equation ( |3.3| ). 

Under these conditions also the parity class of F(to, c, ti, ■ ■ ■ , t^n+i j ) will 

depend only on the parity of c. As in theorem (||) the fundamental polyno- 
mial solution to 

C 2 - F(t ,c, , t Ln ±ij )F 2 = (-If- 1 

is given by C = a Q n + P„, H = Q n . □ 

4. A Worked Example 

As an example, consider those positive integers whose square-roots have 
continued fraction expansion with period of length 9. Thus the symmetric 
part of the period has length 8 and it is necessary to check the 2 4 = 16 
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zero-one sequences to determine which are permissible. (This checking is 
done in essentially the same way as in Example [l] above.) There are 11 valid 
sequences: 

0,0,0,0,0,0,0,0 
0,0,0,1,1,0,0,0 
0,0,1,1,1,1,0,0 
0,1,0,0,0,0,1,0 
0,1,0,1,1,0,1,0 
0,1,1,1,1,1,1,0 
1,0,0,1,1,0,0,1 
1,0,1,0,0,1,0,1 
1,0,1,1,1,1,0,1 
1,1,0,0,0,0,1,1 
1,1,1,0,0,1,1,1 

The ninth of these is considered in more detail (Each of the others can be 
dealt with in a similar way). For clarity the letters a,b,c and d are used 
instead of ti,t2,ts and £4. Evaluating the continued fraction 

(4.1) 

11111111 

° + 2a + 1+ 26 + 2+ 2c + 1+ 2d + 1+ 2d + 1+ 2c + 1+ 26 + 2+ 2a + 1 

it is found that 

P& = Qj = -l-2d + 

2(3 + 4a + 26 + 4a&)(4 + 36 + 4c + 4&c + 6d + 4&d + 8cd + 86cd) 
+ 4(3 + 26 + 4(l + 6)c)(2 + 6 + 3c + 26c + a(3 + 26 + 4(l + 6)c))x 

(l + 2d + 2d 2 ), 

P 7 = 4(l + 6)(4 + 3& + 4c + 46c + 6d + 46d + 8c(i + 86cd) + 

2 (3 + 2 b + 4 (1 + b) c) 2 (1 + 2 d + 2 d 2 ) and 
Q 8 = 8(2 + 6 + 3c + 26c + a(3 + 26 + 4(l + 6) c)) 2 (l + 2d + 2d 2 ) + 
(3 + 4a + 26 + 4a6)(3 + 4c + 4d + 8cd+(2 + 4a)(4 + 36 + 

4c + 46c + 6d + 46d + 8cd + 86cd)). 

Since n is 8 (even) and Q 8 is odd (so to cannot take half-integer values), in 
this case Fg(to, a, b, c, d) is defined by 

(4.2) F 9 (t ,a,b,c,d) = (Q 7 P 7 /2 + t Q 8 ) 2 + 2t P 8 + P? 
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and 

y/F 9 (t ,a,b, c, d) = [Q 7 P 7 /2 + t Q 8 ; 2a + 1, 2b + 2, 2c + 1, 2d + 1, 2d + 1, 



2c + 1, 26 + 2, 2a + 1, 2(Q 7 P 7 /2 + t Q 8 )}, 

this expansion being valid for all a,b,c,d > and all to > —Q 7 P 7 /(2Q 8 ) 
and in particular for all to > 0. In these ranges 

C = (Q 7 P 7 /2 + t Q 8 )Q 8 + Ps, H = Q 8 

gives the fundamental polynomial solution to 

C 2 — F 9 H 2 = -1. 

F 9 (t , a, b, c, d) = (Q 7 P 7 /2 + t Q 8 ) 2 + 2t P 8 + P 7 2 = (1 + tg) mod 4 

so that if (Q 7 P 7 /2 + toQs) 2 + 2io-Ps + -P7 is a square-free number for some 
particular a,b,c,d > and some odd to > —Q 7 P 7 /(2Q 8 ), then 



(Q 7 P 7 /2 + t Qs)Q8 + Ps + J(Q 7 P 7 /2 + t Q 8 ) 2 + 2t P 8 + PfQs 



is the fundamental unit in Q ^ y / (Q 7 P 7 /2 + t Qs) 2 + 2i -Ps + ^7 

5. Mystification, Fermat-Pell polynomials of a single variable 

and more on odd-even 

Clearly it is possible to "mystify" this process by replacing each t{ by some 
polynomial giiti) taking only positive values or by replacing 2ij (recalling 
that the continued fraction expansion contains only terms like 2t{ + 1 or 
2ti + 2) by some polynomial &j(tj) taking only even non-negative values or 
by setting tj = U(Xi, X2, • • • , X^), 1 < % < \ , a polynomial in the Xj's 
taking only positive values, where the -Xj's can be independent variables 
and k can be as large as desired and so on. 

Finally of course one can obtain single-variable Fermat-Pell polynomials 
by replacing the original variables to,ti,l < i < Lhf"J by polynomials in 
a single variable. If it is desired that the period of the continued fraction 
expansion of the new single- variable Fermat-Pell polynomial should stay the 
same as that of the originating multi-variable polynomial then the domain 
of the single variable should be restricted so that the polynomials replacing 
each of the tj's take only positive values as in the multi- variable case and 
the polynomial replacing to must be such that the ao term stays positive for 
all allowed values of the new single variable. 

For example, letting a = s,b = 0,c = s,d = and to = s in the polyno- 



mial (4.2) above produces the single-variable Fermat-Pell polynomial 



g(s) = 639557 + 6858268 s + 33078145 s 2 + 



94534688 s 3 + 177380352 s 4 + 228442240 s 5 + 204593408 s 6 + 

125870080 s 7 + 50925568 s 8 + 12238848 s 9 + 1327104 s 10 
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which has the continued fraction expansion (valid for all s > 0) 



^/gjsj = [799 + 4289 s + 9184 s 2 + 9856 s 3 + 5312 s 4 + 1152 s 5 ; 



2s + 1, 2, 2s + 1, 1, 1, 2s + 1, 2, 2s + 1, 



2(799 + 4289 s + 9184 s 2 + 9856 s 3 + 5312 s 4 + 1152 s 5 )]. 
g(s) = (1 + s 2 ) mod 4 so when s is odd and positive and g(s) is square-free 

51982 + 534625 s + 2429840 s 2 + 6408000 s 3 + 
10812928 s 4 + 12115200 s 5 + 9019392 s 6 + 4304896 s 7 + 1196032 s 8 + 



147456 s 9 + V / g(s)(65 + 320 s + 576 s 2 + 448 s 3 + 128 s 4 ' 

is the fundamental unit in Q[ g(s)]. For example, letting s = 1 gives that 
47020351 + 1537V935888258 is the fundamental unit in Q[ V935888258]. 
Starting with the continued fraction 

11111111 

+ 



2a + 1+ 26 + 2+ c + 2e+ 2d + 1+ 2d + 1+ c + 2e+ 26 + 2+ 2a + 1 
and following the same steps as above with the continued fraction ( |4.1| ) a 
multi-variable Fermat-Pell polynomial is developed which takes values in 
the parity classes associated with permissible sequences 7 and 9. Letting 
a = b = d = e = t = one gets the single-variable Fermat-Pell polynomial 

g(c) = 4325 + 28140 c + 83652 c 2 + 147440 c 3 + 168000 c 4 + 

126528 c 5 + 61504 c 6 + 17664 c 7 + 2304 c 8 
with continued fraction expansion 



y/g(cj = [65 + 214c + 288c 2 + 184c 3 + 48c 4 ; 



1, 2, c, 1, 1, c, 2, 1, 2(65 + 214c + 288c 2 + 184c 3 + 48c 4 )], 

valid for c > 1. 

6. Concluding Remarks 
Every Fermat-Pell polynomial in one variable, s say, that eventually has a 



continued fraction expansion of fixed period length can be found from (3.2), 
if it takes values in only one parity class for all sufficiently large s, and from 
(|3.6[) , if it takes values in two parity class for all sufficiently large s. (Recall 
remark (i) after theorem (||)) 

Of course none of this does anything to answer Schinzel's question of 
whether every Fermat-Pell polynomial in one variable has a continued frac- 
tion expansion. Neither does it provide a criterion (such as Schinzel's in the 
degree-two case) for deciding if a polynomial of arbitrarily high even degree 
is a Fermat-Pell polynomial. Perhaps it raises another question - Does every 
multi- variable Fermat-Pell polynomial have a continued fraction expansion? 
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Does every multi-variable Fermat-Pell polynomial have a continued fraction 
expansion, assuming every Fermat-Pell polynomial in one variable does? 
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